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(Abstract) A statistical study for the switching processes in transmission and distribution electric systems, including earth 
return effect, is presented. The frequency dependence for the parameters of the studied system with a homogenous earth return 
is included as well as the random values for the closing angle of circuit breakers are taken statistically. Expressions for voltage 
and currents in modal coordinates are deduced. The overall presence number of transients and the average over- voltage are 
computed for random angles of switching. The effect of length of lines at different voltage levels is evaluated. A profile for the 
boundaries of either transmission or distribution network is discussed and analyzed. The comparison in the analysis accounts 
the effect of transposition for the non-transposed lines. The results are considered for the standard values of voltage level, which 
is varied from 3 kv up to 500 kv. 
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1 . INTRODUCTION 

Transients are used as disturbance sources for illustration. The 
proposed first method utilises harmonic voltages and wavelet 
coefficients as PQ features. An enhanced genetic algorithm is 
proposed for solving the modified K-means clustering to place 
PQ measurement monitors. This also, would be presented with 
the transients because the assymetry should be presented [1]. 
The transient processes in a power network may lead to damage 
for the operation either partially or completely in spite of it is 
related to the operation system itself. Otherwise, the switching 
transients may be internal or external and consequently, the 
study of such item must be significant to bring the system to the 
safe and suitable zone. Also, the transient phenomena can be 
occurred either for a long time length such as electromechanical 
type or for a short time as the electromagnetic switching 
processes [2]. This subject takes the action with both 
transmission and distribution systems so that it may be 
investigated more and more in future. The impedance Z of a 
symmetrical line as a function of its parameters (conductor 
resistance R c and ground R g as well as self inductance L and 
mutual M may be formulated as given by equation 1 in 
appendix. 



C1 



H1 



D12. 



C2 



D31 



D23 



V 

▲ 



C3 

H3 



H2 



Fig. 1: The geometry of conductors 



There are some different approaches tried to regulate the 
parameters of a power system generally [3-4]. Others went to 
transmission lines or even to change the concept of 
transmission but the final practical solution requires more effort 
[5 -7]. Also, there are new types of lines such as flexible 
alternating current AC transmission systems FACTS [8, 9]. 
The switching transients depends to a great extend on the 
parameters of the system where the ground return effect plays 
an important role. The earth will be considered as a 
homogenous for simplicity and so, the Carson concept would 
be applied [10]. Otherwise, the capacitance is taken as a 
constant without any ground effect while the frequency 

Since shunt capacitors are installed normally in large distribution systems 
and so, many critical nodes may be appeared. Also, the power harmonics and 
capacitor switching 
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dependent parameters are introduced. The self capacitance will 
be C s and the mutual will be C m . Hence, this paper studies both 
systems transmission (500, 330, 220, 110 kV) and distribution 
(35, 10, 6 kV) as sample for any network in order to conclude 
the statistical performance for a system under the transient 
conditions. Hence, the random five values for the switching 
angle are applied to present the general case of switching 
transients [10]. The overall geometry of phases on standard 
towers, as shown in Fig. 1, is listed in Table 1 [1 1]. It should be 
mentioned that the transmission lines of high- voltage networks 
are protected against lightning strikes and internal overvoltages 
with the use of overhead ground wires and surge arresters. The 
appropriate modelling of a power system's components is a 
significant issue, in order to obtain reliable theoretical 
simulation results for the lightning performance of the line that 
will be used in the design of a new line or the improvement of 
an existing one [7, 12, 13]. 

It is required to illustrate the objective aim of the study where 
certain national power system cannot be identical to the other. 
Then, the determination of a general performance for all 
networks through the investigation of any one may take us 
away from the actual condition in the site. However, the 
transient characteristics depend on mainly the parameters of a 
system where their values can be evaluated through the 
geometry of phases. This means that the geometry allocation of 
phases will reflect the overall behavior in the domain of 
electromagnetic transients. Therefore, the classes of 
distribution and transmission can be taken as the gate to solve 
such problem and consequently its generalization may be 
implemented. 

Table 1, a: Dimensions of considered Towers (m) 



Type 


Outdoor 


Indoor 


Tower 


Metal 


Wood 


Metal 




500 


220 


220 




330 


110 


150 


kV 


220 


35 


110 


110 
35 
6-10 


6-10 


35 
20/10 
6/3 




20-19 


12 


5.3 




16-30 


9 


3.5 


HI 


18.5-30 


11 


2.3 


15.5-18 
12-17.5 
9-8.9 


9.5 


0.93 
0.58/0.38 
0.29/0.25 




20-19 


12 


3.5 




16-23 


9 


2.3 


H2 


14-24 


11 


1.5 


12.5-15 
9-14.5 
9-7.6 


8.75 


0.61 
0.38/0.25 
0.19/0.35 




20-19 


12 


1.7 




16-23 


9 


1.1 


H3 


14-18 


11 


0.7 


12.5-12 
9-11.5 
9-7.6 


8 


0.29 
0.18/0.12 
0.09/0.065 



Table 1, b: Dimensions (m) for the considered towers in the 
research 



Type 


Outdoor 


Indoor 


Tower 


Metal 


Wood 


Metal 




500 


220 


220 




330 


HO 


150 


KV 


220 


35 


no 


HO 
35 
6-10 


6-10 


35 
20/10 
6/3 




H-13 


5 


1.8 




8.5-9.4 


4 


1.2 




6.8-7.1 
4.24-3.27 
3-3.08 
1.5 


3 

1.5 


0.8 
0.32 
0.2/0.13 
0.1/0.07 




H-13 


5 


1.8 




8.-8.6 


4 


1.2 


D23 


6-6.86 


3 


0.8 


5.5-3.27 
3-3.06 
1.5 


1.5 


0.32 
0.2/0.13 
0.1/0.07 




22-26 


10 


3.6 




16-7.5 


8 


2.4 


D3l 


6.8-13 


6 


1.6 


4.24-6 
4.24-6 
3-1.5 


1.5 


0.64 
0.4/0.26 
0.2/0.14 



This projects the idea of the presented work according to the 
evaluation of the boundaries for the profile of switching 
transients in power system. The investigation may be tailored 
into two steps. The first will be for each system (Distribution or 
transmission) while the second deals with the statistical 
variation in length or voltage or even in the percentage content 
of a system relative to the other inside the network. After that, 
the general profile can be computed for different values of 
content ration (Distribution to transmission ratio). Otherwise, 
the data cannot be based on the population readings so that the 
statistical sampling processes must be inserted. 

The main equations of voltage and current at a point x of a 
line are given in a matrix form of equation 2 in the appendix. 
Also, the deduced Voltage at a point x with the help of 
Sylvester theorem [12] may be found as shown in appendix 
(equation 3). Otherwise, the current can be derived as in 
equation 4 (Appendix). The replacement of the matrix of 
propagation coefficient by a multi term equation according to 
the rule of Sylvester will simplify the derived formula for the 
propagation coefficient (See equation 5 in the appendix. Then, 
the formula of the propagation coefficient appears to be in the 
form of equation 6 (Appendix). 

2 . MODAL COORDINATES {a, p, 0} 

This result will be more complicated for the non-transposed 
lines in the LAPLACIAN domain and so the final formula can 
be solved in the wave mode coordinates. On the other hand, the 
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distribution power networks plays an important role in the 
processes of operation and many papers were analyzed this 
system as in [4]. The modal technique gives the simplicity for 
the solution of such problems in the complex plane so that the 
main parameters of the deduced equation (6) may be applied in 
each of the modes of the analysis. 




SDW 3BDW 23W . TKW 35W 101V 

Fig 2 The computed values for the capatitance in the wave 
modes of different lines 

It should be mentioned that, the wave modes are individual 
channels without any mutual relationship between them and 
they are defined as {(a), (P), (0)} as deduced in appendix [10]. 
It is marked that the speediest wave mode is the first (a) while 
the second (P) is slower a little. The last important wave mode 
is the zero wave mode (0) which is depending on the earth 
return effect. This earth is normally a non homogeneous 
medium either vertically or horizontally. This item has a great 
interest of the research field because all unexplained 
phenomena may relate to it. It is seen that, both first and second 
wave modes in modal coordinates are practically the same and 
varying in the same manner although the zero mode value can 
be considered a constant for all levels of voltage. 

Therefore, the parameters of the system must be evaluated. 
So, a capacitance is calculated for different phase geometry 
shapes in modal coordinates as seen in Fig. 2. Its value is 
increased in the second wave mode for the Extra High voltage 
(EHV) and ultra High voltage (UHV) while it still a constant 
for the High voltage (HV) level till the distribution system. In 
the distribution system the effect of the capacitive reactance on 
the presence of transients is decayed after the 100 kV level and 
before the 35 kV level. It begins the decrease process before the 
35 kV level. 

It is important to indicate that the transient analysis, 
generally, contains both time and frequency domain where they 
sometimes are mixed together or otherwise, the transformation 



from one to another may be required. In some methods the 
procedure of solution tends to model a three-phase transmission 
line directly in time domain, without the explicit use of inverse 
transforms for simplicity. This must take into account the 
frequency-dependent parameters of the line, considering both 
the soil and the skin effects. Furthermore the accuracy of the 
developed results should be verified, in frequency domain, by a 
simple methodology based on line distributed parameters and 
transfer function related to the input / output signals of the 
lumped parameters representation. The commercial software of 
electromagnetic transients program (EMTP) may be 
considered for check of results of new methods. 
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Fig 3 Variation of propagation coefficient for the 
studied lines(for voltage level of 500&220kv) 

Two limits (upper and lower) are indicated where both 
channels { 1 st (a) & 2 nd (P)} have a higher value than that for the 
zero wave mode (0). It is increasing with voltage level in the 
two modes but it is a constant in the zero wave mode (0). The 
propagation coefficient is varied in these coordinates a shown 
in Fig. 3 where per unit system is considered for the value. Also, 
the variation has been computed according to the frequency in 
the frequency domain. 

The frequency presence inside the processes of transients 
may reach 60 kHz due to the high frequency HF loads and 
interference with electromagnetic fields. The 50 Hz t zero mode 
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position is taken ass the reference (unity value) and then its 
value is illustrated for studied voltage levels. The coefficient is 
increasing with frequency in Zero mode coordinates while it is 
decreased in other channels. The maximum ratio of this 
coefficient is appeared for 1 10 kV, where the 35 kV distribution 
lines approaches the case of 500 kV. This approves that the 
distribution network is subjected normally to heavy transients 
relative to its level where its value may exceed that for EHV 
levels. This coefficient can be shown mathematically in each 
channel from equation (6) to obey the formula 7 in the 
appendix. 

On the other hand, the propagation factor, as shown in Fig. 3. 
illustrates that there is a difference between both first and 
second wave modes where the first wave mode has a maximum 
value for all levels of applied voltage. This indicates that the 
first wave mode is the speediest one while the zero wave mode 
has the slowest. Contrary, the zero mode parameter begins at 
lower value which increased with the rise of frequency. The 
frequency means the presence of harmonics in the real wave 
forms. These forms can not be sinusoidal as the normal 
operation conditions but the increase of higher frequency 
means more harmonics. The propagation factor becomes 
constant quickly and becomes frequency independent in all 
three modes. 

However, the attenuation factor plays great role in the 
limitation of the over-voltages since it is representing the 
effective part of the propagation process. The characteristic of 
the average transient voltage calculated with respect to the 
nominal operating voltage is shown in Fig. 4 where the 
maximum ratio is corresponding to the higher voltage class. 
The PER UNIT system is considered where the nominal 
operating voltage has been the reference value for each one. 
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Fig 3 (b) Variation of propagation coefficient for the 
studied lines(for voltages level of HOkv) 
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Fig 3 (c) Variation of propagation coefficient for the 
studied lines (for voltages level of 35kv) 



It should be indicated that, the attenuation factor is very high 
for UHV and EHV but is highly decreased in the distribution 
systems after the 35 kV level networks. Also, the HV level 
gives a constant effect for the value of attenuation factor where 
this region begins at 220 kV up to the distribution level of 35 
kV. 

However, the wood tower standards give a lower attenuation 
so that its value for 35 kV towers may exceed that for the 110 
and 220 such towers. The derived equations will be more 
complicated for the non-transposed lines in the LAPLACIAN 
domain (p) and so the initial formula can be solved in modal 
coordinates. Also, the distribution networks play an important 
role in the operation conditions s analyzed before [8]. Many 
methods are known for the transient calculations and the wave 
mode propagation method will be selected because it is valid 
for either transposed or non-transposed lines. It mainly depends 
on the transformation matrix [T], which takes the form of 
equation 8 in the appendix. 

If the basic differential equations (2) are specified to the 
modal coordinates, new deduced expressions become 
individually (equation 9 of the appendix). Then, voltage and 
current at point x in modal coordinates will be in the form of 10 
in the appendix. These equations depend on the transformation 
of the phase system into another one, which is known as the 
wave mode system. This means that the wave modal concept 
has three isolated modes without any mutual effect between 
them as this can be reached through the mathematical theorem 
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of Eigen values and Eigen vectors. It should be noted that the 
transformation matrix of equation 8 is the Eigen vectors for the 
system. 

A special program is used for the calculations of parameters 
in wave modes or the transients and modified to measure the 
statistical values as presented here. The application of 
parameters in modal axes is based on the matrix equation 11 
(See Appendix). 




Fig 4 (a) The mode performance for average transient 
voltage at different nominal operating voltages 

For easy computation these equations must be referred to the 
voltage source at the sending end instead of the voltage t 
receiving end. Then, a simplified formula for voltage and 
current (equation 12 of Appendix) can be derived (See 
Appendix). It would be indicated that the voltage source differs 
from that at the sending end of line according to the relation 
(equation 13 in the appendix). Consequently, the voltage may 
be formulated through the transformation matrix (See formula 
14 in the appendix). Therefore, voltages in phase coordinates 
can be determined numerically according to the convolution 
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theorem [12] as given in appendix (equation 15). Similarly, the 
equations of currents in phase coordinates can be found [12]. 

It is seen from Fig. 4 that, the maximum value for the 
standard deviation is appeared at the distribution systems where 
it comes at the voltage level of 35 kV. This indicates that the 
distribution level may be also, affected by the switching 
processes. This leads to the importance of the use of the 
capacitive banks in the distribution level. So, the automated 
units for compensation in the distribution systems would be 
checked for all conditions because the deduced results give a 
danger pointer for the presence of transients in such systems. 




Fig 4 (b) The calculated standard deviation at different 
nominal operating voltages 

3 . STATISTICAL SWITCHING TRANSIENTS 

Moreover, the variation of time domain length has been 
examined as given in Fig. 5. Interval duration of 0.21 ms was 
considered for 100, 200 and 400 points of calculation. The 
attenuation effect is drawn in Fig. 5 in frequency domain while 
Fig. 5, (a) illustrates the variation in the first mode. Fig. 5, (b) 
shows this dependency for the most important wave mode 
which is the zero wave mode. Its importance comes from the 
slow of traveling waves according to the natural values of the 
elements and components in the transmission lines with 
distributed parameters. It should be mentioned that, the 
consideration of asymmetry not only for the distributed 
parameters but also the absent of transposition between phases. 

It is seen from Fig. 5 that, the distribution system at the level 
of 35 kV has a great values that exceed the transmission system. 
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Fig 5 (a) average over- voltage at the receiving end of lines for 
different values of intervals 

calculations are based on the voltage at the receiving end of 
transmission lines. The results prove that the small interval for 
the calculations gives a rise in the average transient voltage and 
consequentially, the evaluated number of transient voltage. 

In the above results, the combination of changing of either 
interval or number of points in the domain length clarified the 
effect of steady state presence in long time term for 500 kV. 
The low voltage distribution system (6-10 kV) gives gradually 
decrease with time domain length while it is approximately 
constant for 35 kV. Thus the choice of the small interval of 
calculations must be considered, and accordingly a suitable 
duration for the overall time of calculations should be taken. 
This means that the steady state time after switching processes 
would be eliminated as possible in order to get real exact 
results. 
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Fig 5 (b) attenuation factor for waves at the receiving end 
of lines for different values of intervals (zero mode) 

Also, the overall average voltage for studied lines is 
decreased due to the presence of steady state period inside the 
process of calculations. Otherwise, the effect of interval only 
for a constant value for the points of computations (300) is 
considered and the results are listed in Table 3 for both verge 
voltage and standard deviation. 




Fig 6 (a) Number of over- voltages at the receiving end 
of lines for different values of intervals (%) 

Therefore as it is concluded that, the choice of time interval 
for the computational analysis affect on the results of proposed 
number of transients so that an optimal value for the chosen 
interval would be pre- selected before the study of any condition. 
Also, this number will be a good measurement for the presence 
of over- voltages inside a process where the sever cases must be 
avoided. The planning for the transmission systems as well as 
the distribution systems should account the phenomena of 
transients so that its consideration may be a vital parameter for 
the design body. The transmission systems take the major 
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importance in this presented occasion because the transient 
level will be very high so that its effect may be accounted as a 
first step in the planning scheme. This effect may be 
disappeared in the distribution level at the consumer use. 

1.4 j- 
1.3 
1.2 



1.1 



21 ms32.2 ms42 ms 69 ms 84 ms 

Fig 6 (b) Number of over- voltages at the receiving end of 
lines for different values of intervals (%) 

4 . TRANSIENT NUMBER 

The transient presence in a system during the switching 
processes may be necessary to measure the level of its action on 
the insulation level. This paper presents a new proposed 
number for such simulations in order to measure the transient 
content TN in a system. It represents the percentage presence of 
over- voltages inside the time length considered and so, it can be 
expressed mathematically in the form: 



TN = N t / K 



(1) 



The presence of transients for each case is drawn in Fig. 7, a, 
but the overall presence will be shown in Fig. 7, b. It is seen that, 
the transient level is always decreasing with voltage level 
except some readings of 220 kV relative to 110 kV due to the 
closed dimensions of towers for both levels (Table 1). This 
voltage is decreased with respect to the interval, which reflects 
the total time domain length. This confirms the deduced results 
above of Fig. 6. Also, the same conclusion can be derived from 
Fig. 7, b. 
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Fig 6 (c) Number of over- voltages at the receiving end of 
lines for different values of intervals (%) 
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Table 3, a: Average voltage level and standard deviation a at 
receiving end R E of lines for different duration of intervals at 
K=300 (ms) 



kV 


23 (ms) 


22 (ms) 


V 


a 


V 


a 


500 


1.221 


.007315 


1.23 


.007373 


330 


1.129 


.00846 


1.135 


.00902 


220 


1.069 


.00739 


1.075 


.007491 


110 


1.072 


.01488 


1.085 


.0147 


35 


1.075 


.02439 


1.065 


.02348 


6-10 


1.014 


.006732 


1.015 


.007079 



Table 3, b: Average voltage and standard deviation a at the 
receiving end of lines for duration of intervals at K=300 (ms) 



kV 


21 (ms) 


12 (ms) 


V 


a 


V 


a 


500 


1.23 


.007972 


1.291 


.008629 


330 


1.135 


.008887 


1.218 


.009699 


220 


1.075 


.008276 


1.123 


.00754 


110 


1.085 


.0162 


1.081 


.009005 


35 


1.065 


.02144 






6-10 


1.015 


.006912 







This table lists the average voltage as well as the standard 
deviation for each studied voltage level with respect to the 
duration of calculations. So, both standard deviation value and 
the average transient voltage are decreased with the increase of 
the duration of calculations in spite of the small difference 
between all durations. 

Fig. 7 ensures the above conclusion corresponding to the 
long time choice where the steady state condition may be 
reached. This is appeared at the time of 63 ms the transient 
contents become less than 1.1 although it is very large value for 
36 - 60 ms instead. This value will be decreased a little in the 
range of 60 - 63 period of the calculation. Thus, a small time 
not exceeding the calculation period of 60 ms may be proposed 
for applications. 
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Fig 7 (a) The result for the studied lines at different 
duration intervals at(K=300) 
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Fig 7 (b) The result for the considered lines at different 
duration intervals at (K=300) 

In other words, Fig. 7 presents a large value of percentage 
content of transients in the system is appeared at the 500 kV for 
all dimensions of this class of voltage where it varies between 
1.35 and 1.2. This value is steeply decreased for the distribution 
system where the transient switching voltages are computed as 
1.01 up to 1.03. In a midpoint lies the transmission of 330 - 1 10 
kV level where its content varies between 1.2 and 1.07. 



5 . COMPENSATION EFFECT 



The reactive part of elements connected to a power system 
takes a part in the switching processes as soon as voltage 
control in the distribution system needs shunt capacitor on 



loads. This may be transferred (directly or non-directly) to the 
low voltage line in the case of heavy loads. This capacitor may 
be appeared also in EHV and UHV transmission systems. 

The design of transmission systems depends to a great extent 
on the computational analysis of switching transients because 
the generated voltages in switching processes will be danger so 
that the insulation level must be stated to be able to withstand 
with these high voltages. Contrary, in the distribution systems 
the switching processes may be at a less level of importance 
although this is not the general rule as it is appeared in the 
above results. This may be seen from the value of the calculated 
value of over-voltage in the system where it is generally 
permits the use of capacitor banks for the power factor 
corrections. 

Whatever, the load type has been tested where both inductive 
and capacitive loads are introduced at the receiving end of the 
switched lines. Results are shown in Fig. 8 for the inductive 
condition for a time length of 200 points with .023 ms interval. 
It shows a normal characteristic while the heavy transients 
appear with the second condition (capacitance). 
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Fig 8 (a) The average value of voltage for the inductive 
load effect on transient processes 



It is very important to avoid the presence of the capacitive 
reactance in the elements of the power system because it has a 
danger effect on the recovery voltage. This may be expressed 
here through the presence of transient voltages issn the system 
where it may be reflected by the transient number. 
Consequentially, the harmonics may be too generated due to 
the presence of capacitance in the system. So, the nonlinearity 
characteristics could become more complicated in the solution 
in closed form. 
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Fig 8 (b) Standard deviation for inductive loads 

It would be required to indicate that the capacitive effect 
reached faster in the distribution system. It began at 10 Ohms 
for 6 kV with a per unit value of 37. 14 wshile this value became 
16.82 for 35 kV lines. The other lines still in the normal but the 
resonance condition may come later. It varied with potential 
level as the final will be the highest one. 
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Fig 9 (a) Average Transient voltage for capacitive loads 

The effect of capacitive load on the transient profile could be 
a danger because it contains a lot of harmonic of high orders. 
This means that resonance characteristics may be appeared due 
to switching processes 

The resonance characteristics are shown through curves in 
Fig. 9 although the standard deviation became very large. This 
means that the value of actual maximum voltage is high. The 
curves for 220 and 330 kV are not illustrated. 
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Fig 9 (b) Standard deviation for capacitive loads 



Fig. 10 shows well that the presence of transient inside the 
time length is higher with nominal potential and it is very high 
for the capacitive condition. The average voltage distribution at 
the receiving end of the lines studied but for all possible lengths 
is computed as given in Fig. 11. 
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Fig 10 The overall transient number in both cases of 
inductance and capacitance 

It proves that the switching transients in EHV and HV lines 
will have a great role in the design of such lines and it decayed 
for the distribution system. Also, the effect of line length shows 
an average value for each length in spite of the different 
potentials (Fig. 11, b). 
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Fig 11 (a) The calculated average value for different line 
lengths at K=200 



The derived standard deviation is drawn in Fig. 12 with 
different lengths of lines when the transient number is 
presented in Fig. 13. These curves indicated the same 
conclusion defined above. 
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Fig 11 (b) The calculated average value for different 
line lengths at K=200 



Nevertheless, the shunt compensation by the installation of 
rector at receiving end would be the best solution for long 
lines EHV and UHV systems due to the reduced statistical 
transient presence (as the value of TN). Whatever, shunt 
capacitors on feeders of thee distribution networks 
represent critical condition with bulk reactive powers so 
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that a great interest would be needed with sudden change of 
load according to the performance of load curves at such 
points. 

6. STATISTICAL PROFILE 

The process of profiling for a group of data may be expressed 
mathematically either through a deterministic or by a stochastic 
concept. The deterministic style is a mathematical model, 
which means the actual specified all data are included in the 
system while the stochastic concept gives an overall view for 
the system through either all data or some random values. 
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Fig 12 The calculated average value for different line 
lengths at K=200 

Although this last concept may introduce an error according 
to the statistical principles, a good result expressing a problem 
can be easy achieved [13, 14]. This will be more difficult with 
the internal transients. Therefore, the above analysis can be 
used to get a general profile for the processes of switching 
transients in a certain specified network. 

Any network contains different numbers of distribution and 
transmission systems with various connections [15, 16]. There 
is a specified ratio of these systems to a certain moment but this 
ratio will be varied with each extension, which must be 
normally happened, continuously, at least in the distribution 
system. Nowadays, new concepts of modern techniques would 
be needed so that the Artificial Neural Network (ANN) can be 
implemented for the modification of the operation in electric 
power networks [17 - 19]. So, the ratio between both systems 
inside the network will be actually a dynamic ratio. Therefore, a 
general profile for such specified network at a certain moment 
can solve the problem of generalization of the view for this 
network. Thus, the distribution to transmission ratio (DTR) my 
be defined by 
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Fig 13 The deviation transient for different line lengths at 
K=200 



(DTR) = E(Z,V Z(Z,) r 



(2) 



Consequentially, line impedance can be computed according 
to the specific resistance and inductance multiplied by the 
length and then, the summation of impedance for either 
transmission or distribution would be evaluated. After that, a 
defined profile for a network should be deduced for at a certain 
network at a specified time but a general profile for this 
network may be evaluated at all values of (DTR). 

On the other hand, a general ratio (GR) may be proposed as: 
(GR) = DTR / (1+ DTR) (3) 

Table 2: The Statistical proposed value for the percentage 
content of transients in: 
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Fig 14 The deduced profile for the proposed cases of 
transmission and distribution system 

This will be useful to find the general profile for a network 
through the general transient number by: 



TN G = (GR)[ (N t ) D -(N t ) T ]/K + (N t ) T /K 



(4) 



This general transient number may be evaluated for system 
specified, but it may be important to find its change due to the 
growing in the electric systems as a whole. So, a basic 
configuration for a power system is suggested as listed in Table 
2 where both transmission and distribution systems re indicated 
according to the self parameter for each voltage level. 

Then, different cases with the variation of such power system 
are considered on the bases of normal growth in the developing 
counties. The impedance is referred to the 1 1 kV impedance as 
calculated statistically above while the percentage content for 
either the transmission (500, 220, 110 kV) or the distribution 
(35, 11, 6 kV) level with the determined transient number. 
Therefore, the constant configuration for the transmission 
system is considered (basic case), where the variation is applied 
to the distribution level. The general transient number is 
computed for each case as shown in Fig 14. 

Also, the change in the transmission system is given (cases 4 
and 5). It is shown that the transient content in a network will be 
increased with the growing in the distribution network. This 
may give more increase with the great extension in the lower 
level of voltage so that a corresponding extension will be 
needed as proved for the case 4, where the 1 10 kV increased. It 
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should be indicated that this extension balance cannot be 
achieved t the highest UHV scale (case 5). 
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Fig 15 The determined rate of rise of the transient content 
referring to the base case for the proposed cased of 
transmission and distribution system (%) 

Thus, the proposed technique will be useful for the planning 
and design of power systems in order to decrease the transient 
level inside a network [20]. This may be clarified with the rate 
of rise of the transient content as given in Fig. 15 for the studied 
cases referring to the basic one. The rate of rise is increased 
with 6 kV (35 %) although it reaches 40 % for the UHV 
extension. This proves that the proposed concept for the general 
transient number must be considered for any extension to avoid 
the heavy presence of transients in a power system. 

7 . CONCLUSION 

This work reported a statistical evaluation of the switching 
processes in transmission and distribution systems with earth 
return effect. The frequency dependence of the parameters was 
investigated and the random values for the closing angle of 
circuit breakers were statistically taken as well. 

The paper has proposed a simplified concept for the 
evaluation of transients in both transmission and distribution 
systems through the transformation into the so called the wave 
mode coordinates. Also, it is valid for both transposed and 
non-transposed phases. 

The use of bulk capacitors in the distribution power network 
has been tested for the transient presence in a percentage style 
with the limits of resonance phenomena. 

The overall transient number for the power network in a 
statistical base has been realized so that the UHV and EHV still 
have the heavier transient presence. 

The statistical profile is a good tool for the planning and 
design of power systems especially with the quick developing 
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networks. This case may be appeared in the developing 
countries where the great density of population will be 
presented. 

APPENDIX 

The presented method here depends on the transformation of 
the problem in the time domain (with mutual effect between 
phases) into a new frequency domain (3 independent modes). 
Then, the solution will be easy to done where the final solution 
in frequency domain must be returned to the time domain 
expression using the inverse transform technique. This process 
will be based on the deduced general transformation matrix [T]. 
It is valid for all types of lines either transposed or 
un-transposed transmission lines. 

Moreover, the asymmetry of the problem either in 
transposition or the return path can be considered. This return 
path is the earth return effect where it is simulated 
mathematically and introduced in the values of the parameters 
such as: resistance R, capacitance C, potential coefficient, 
inductance L , surge impedance Z c , attenuation factor a, and 
propagation coefficient y. These all parameters are investigated 
in the present paper above. A differential equation in time 
domain becomes three linear simple equations in the wave 
mode coordinates (a, (3, 0) where their inverse will be simple as 
given in this Appendix. 

The impedance Z of an electric system in its real R and 
imaginary X parts may be expressed in the matrix form by 
[Z] = [R]+j[X] = 



R c +Rg 

R c 



The return effect is introduced through the ground earth 
resistance where the Carson analysis would be considered. It 
assumes that the ground resistance is a frequency dependent 
parameter. 

Thus, the general differential equation of a transmission line 
could be stated as 
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-d/dx [V (x)] = [Z] [I (x)] 
-d/dx [I (x)] = [Y] [V (x)] 



(2 a) 



The solution of these equations can be deduced for the voltage 
V(x) at a point x on the transmission line in the form: 

V (x) = (l/3){ch Yl (l - x) [Mj] + ch y 2 (l - x)[M 2 ]) [v (1)]}+ 
(1/3){Z! sh y 2 (1 - x) [M 2 ] + Z 2 sh Vl (1 -x) [MJ) [I (1)]} (3 a) 

Where the matrices of constants Mi and Mi are formulated as 



[M 1 ] = 



11 1 
11 1 
11 1 



[MJ = 



2 
1 
VI 



-1 -1 
2 -1 
-7 2 



Also, the impedances Zi and Z 2 can be expressed by 
ZiMR + ptL-JVDI/lp (C- C m )} 



Z 2 2 = {R + 3R + p (L + 2M)} / {p (C + 2 C m )} 



Similarly, the current I(x) at a point x would be deduced by the 
expression 

[I (x)] = {sh y 2 ( / - x) [ZJ 1 [v (I)] + ch Yl ( / - x) [I (I)] (4 a) 

The propagation coefficient y can be formulated in the matrix 
form 



[y] = (l/3)Ri [M 1 ] + X 2 [M 2 ]} 



(5) 



It should be mentioned that the value of X variables would be 
given in the relationships: 



X l 2 = y 1 + 2y 2 &X 2 2 =X 3 2 



The expansion of these relationships can be estimated in the 
form: 



e [Y]x = (l/3){e\ x [Mi] + e V [M 2 ] } 



(6) 



e M x = *{(y p - y) (y - y)}/{(yp - y a ) (y - y a )} +e y P x {(y a - y) (y ■ 
Y)} / {(Ya - YpXYo - Yp))+e Y o* {(Ya - Y)(Yp - ?)}/{( Ya - Yo) * (Yp - 
Yo)}= 
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For the non- transposed line^J the transformation matrix is 
expressed in the general form as ~~ 



[T] = 



This matrix [T] is the transformation matrix from the time 
domain analysis into the frequency domain expressions 
although its inverse [T]" 1 may be needed for the inverse 
transform in the opposite direction of transformation. It is the 
transformation matrix for the modal analysis. Its special case 
for the symmetrical lines gives the well known matrix "Clarke 
Matrix" where A = 1 and B = 2. 

Thus, the equations should be solved in the wave modal 
coordinates (a , P , 0) as given below: 



-d/dx [V (x) a ,p, ] = [Z] [I (X) a , p , ] 

-d/dx [I (x) a ,p, ] = [Y] [V (x) a , p , ] 
These equations may be simplified as 



(9) 



[V(x,p)] a ,p, =[chy (l-x)] a ,p, [V(l,p)] a ,p, + [shy (1 - 

x )] a , p , [ Z c ] a , p , [I (1? P)l a , p , 

[I(x, p)] a , p , o = [sh y (1-x)] a , p , o [Z c ] _1 a , p , o [V G, p)] a , p , o + 
[chyG-x)] a ,p, [IG,p)]a,p,o (10) 



In the above expressions the characteristic impedance Z c is 
appeared in the wave mode coordinates. The inverse of the 
transformation matrix may be appeared with the impedances in 
the form 

[E] = [T] [Z c ch y 1 / sh y 1] [T] _1 [I (0)] (11) 

So, we can get the voltage V(x,p) and current I(x,p) in the 
frequency domain in a closed form as 

[ V(x, p)] a , p , = [ch y x] a , p , o [ V(0, p)] a , p , o+[sh y x] a , p , 
[ Z c ] a ,p, [I (0, p)] a ,p,o 

[I(x, p)] a , p , o = [sh y x] a , p , o [Z c ] _1 a , p , o [V(0, p)] a , p , + [ch y 
A J a , p , [i (0, p)] a,p,0 

(12) 

Both the supply voltage E (p) as well as the supply current 
Lsuppiy must be introduced in the mathematical analysis using 
the equation 

[V (0, p)] «, p , = [E (p)] „, p , o - p [L supply ] (13) 

Consequentially, the voltage V(x, p) at a point x of the line may 
be expressed in wave mode coordinates [a , p , 0] by 

[V(x,p)] „ , p , o = {[ch y (l-x) /sh y Z] „ , „ , [E(p)]} /{B 

[ Z c ] a {[ch y (l-x)/sh y Z] „ +A [ZJ P {[ch y (l-x) /sh y Z] „ } (14) 

Finally, the voltage of phases V a ,b, c ( x > P) in phase coordinates 
(a, b, c) will be derived as 

V a (x, p) = AV a + BV p +(A + B) V 
V b (x,p) = AB(V p -V a ) 

V c (x,p) = AV a + BVp-(A + B)V (15) 
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